Is relaxation to equilibrium hindered by transient dissipative structures in closed 

systems? 
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Dissipative structures are generally observed when a system relaxes from a far from equilibrium 
state. To address the reverse question given by the title, we investigate the relaxation process in 
a closed chemical reaction-diffusion system which can potentially form Turing-like patterns during 
the transient. We find that when certain conditions are fulfilled the relaxation process is indeed 
drastically hindered, once the pattern is formed. This slowing-down is shown to be due to stepwise 
relaxation, where each plateau in the relaxation process corresponds to residence at a certain spatial 
pattern. Universality of the phenomena as well as their biological relevance is briefly discussed. 

PACS numbers: 



In open systems far from equilibrium, organized struc- 
tures are well-known phenomena^]. Such "dissipative 
structures" include temporal rhythms and spatial pat- 
terns in chemical reaction-diffusion systems, hydrody- 
namical systems, optical systems and so forth[ulalalll||. 
Among others, biological complex structures can also be 
maintained by non-equilibrium conditions^], 0. 

In the study of dissipative structures, systems are gen- 
erally prepared in far from equilibrium states by imposing 
certain constraints. For example, the concentrations of 
some chemicals are fixed at (compared to the equilibrium 
state) higher or lower levels by supplying or removing 
them from the outside. On the other hand, in biological 
systems, non-equilibrium conditions are maintained au- 
tonomously, at least when considering long time spans. 
As a first step for understanding the autonomous sus- 
tainmcnt of biological non-equilibrium conditions, it is of 
interest to investigate the possibility that the longevity of 
the conditions for the dissipative structures is extended 
by the formation of the structures themselves. 

In closed systems, of course, equilibrium states with- 
out any structures are reached eventually. However, os- 
cillatory behaviors or spatial pattern formations can be 
observed as transient phenomena durin g th e course of re- 
laxation to equilibrium [1 Fl IsL H Hfl 111 Til^ . Here, we ad- 
dress the following question: Can the formation of (tran- 
sient) dissipative structures make far-from-equilibrium 
conditions last significantly longer by slowing down the 
relaxation process to equilibrium? To answer this ques- 
tion, we study the relaxation behaviors of a closed cou- 
pled chemical reactor that can potentially form transient 
Turing- like patterns during the relaxation process. 

Here, in contrast to most studies in reaction-diffusion 
systems, we need to take the changes in the concentra- 
tions of all chemicals into account, instead of keeping 
the concentrations of some chemicals constant. Thus we 
consider the following reaction-diffusion system consist- 
ing of the reactions I) A + v + 2u t^^a B + 3u ' n ) u ^ v > 
and lll)A + u ^ A + v, and also diffusion. Considering 



the limiting case kAB = 1 >> ksA, the evolution of the 
chemical concentrations is given by 

A,v,uJ 



u, 



(1 + Ai)(ui - Vi), 



(1) 



-AiViuf - (1 + Ai)(vi - Ui) + D v (v l+ i + Vi-i - 2vi), 

(2) 



At 



-AiViuf + D A {A l+1 + - 2Ai). 



(3) 



Here, Ui, Vi, and Ai denote the concentrations of the 
chemical components and i denotes the index of each 
site in a one-dimensional space, where periodic boundary 
conditions are adopted for i. Each chemical diffuses to 
neighboring sites with a diffusion coefficient Dx (X = 
u, v, or A). Although we adopt a spatially discrete system 
for simplicity, the conclusions drawn do not change even 
when the continuum-limit (partial differential equation) 
is taken. The diffusion coefficient D u is assumed to be 
slow, and we mostly study the case with D u = as in 
eq(l) since this will not affect our findings qualitatively 
as long as D u « D v . The chemicals Ui, Vi and Ai 
can be regarded as activators, inhibitors and resources 
respectively. 

This model is a variant of the Gray-Scott 0, E3 or 
Brusselator^l ^| models such that we additionally in- 
clude changes in the resources Ai. Note that the value 
2^ Yji ( u i + v i) — ^ is conserved due to the system being 
closed (N is the number of sites.). 

While relaxation to a unique equilibrium state satisfy- 
ing ^4 = and Ui = i>; = S is assured for t — > oo in this 
model, if we fix Aj = Ao >> 1 in order to maintain the 
non-equilibrium condition, this system shows the follow- 
ing bifurcation of the attractor, depending on S. I) If 
S < 0.75, a unique uniform state with Ui and Vi constant 
over i and time exists that is stable against small pertur- 
bations. II) If S > 0.75, the uniform state is unstable, 
and the attractor is replaced by a non-uniform pattern of 
Ui and Vi, which is constant in time. This Turing insta- 
bility of the uniform states is straightforwardly obtained 
by linear stability analysis. 
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In order to study the relaxation process from a non- 
equilibrium state to the homogeneous equilibrium state, 
we now investigate the effects of changing A^ as in eq.(3) 
for the case S > 0.75. Here, depending on the initial 
configurations of Ui, spatial patterns can be formed dur- 
ing relaxation to the homogeneous equilibrium state. We 
study typical relaxation behaviors by varying initial con- 
figuration of Ui = with S — 4 and the initial condition 
M = Mm = 100 and Vi = S = 4. We control initial spa- 
tial inhomogeneity of by taking an initial condition 
S + S x rndi, with rnd-i as a uniform random number 
over [-1,1] (0 < 5 < S). 

Two sets of typical temporal evolutions of ui and Ai 
are displayed in Figs. 1(a) and 1(b) where S = 0.1 in 
(a) and 5 = 4.0 in (b) with D v = 250 and D A = 0. 
The pattern is plotted until it has nearly reached the 
equilibrium state. The corresponding time evolution of 

< A > (= jf Ef Ai) is plotted in Fig. 1(c). 

When 5 is small, Ui remains almost flat with only mi- 
nor fluctuations, and no structure is formed as in Fig. 
1(a). In this case < A > decreases smoothly with time 
as indicated by the solid curve in Fig. 1(c). On the other 
hand, when 5 is large, the initial inhomogeneity in Ui is 
amplified leading to the formation of a nonuniform spa- 
tial pattern that is sustained over some time span until it 
is re-organized into a different pattern, as shown in Figs. 
1(b) and (d). In this case, the relaxation of < A > ex- 
hibits some plateaus as shown in Fig. 1(c), and requires 
much more time as compared to the case when S is small. 
In general, several plateaus are observed during the re- 
laxation process, each of which corresponds to a specific 
spatial pattern, as shown in Fig. 1(d). 

Hence, we have found an explicit example in which 
the formation of a dissipative structure slows down the 
relaxation process. This behavior is rather general in our 
model, as long as S and D v are large enough to allow for 
the formation of spatial patterns. 

In order to obtain insight into the relationship between 
pattern and relaxation, we have measured the spatial in- 
homogeneity of Ui defined by F(t) = ^ \ u i+i ~ u i\ 2 - 
In Figure 2 we plot the decay rate of < A > defined by 

< A >'= dl ° 9 f t A> , as a function of F{t). As can be seen, 
the system alternates between structure formation where 
Ai is consumed and residence at the formed non-uniform 
structure where consumption of Ai is suppressed. Indeed, 
the decrease of < A >' is highly correlated with the in- 
crease of F. Thus, the slowing down of the relaxation 
process by the spatial structure is confirmed. 

Next, we study the conditions for this slowing down of 
the relaxation process. We investigate the dependence of 
the relaxation time on the initial inhomogeneity. Figure 
3(a) shows the sample average of the relaxation time T as 
a function of the initial heterogeneity S, computed up to 
the time when < A > has decreased to 0.1 Ai n i. Here, the 
parameters are set to D v — 250, S — 4, Ai n i — 100 and 
N = 200, while the diffusion constant Da is chosen to be 



0, 0.25, and 25. As can be seen, there is a critical inho- 
mogeneity S c (« 0.5), beyond which the relaxation time 
increases, when Da is small. Indeed, S c is nothing but 
a threshold for the inhomogeneity above which the the 
re-organization of the spatial structure is possible. For 
large Da, however, the re-organization of the structure is 
even then not possible, and the relaxation time is insen- 
sitive to the initial fluctuations. The threshold S c exists 

for Da ~ 20, while the value of S c itself is insensitive to 
the value of Da within this range. 

The dependence of the relaxation time on Da, on the 
other hand, is plotted in Figure 3(b) where S = 0.1 or 
= 4. When S is smaller than <5 C , the relaxation time 
remains short. In contrast, the relaxation time shows a 
peak around Da ~ 0.3, when S is larger than <5 C , while 
for smaller and larger Da, it approaches constant values. 
Here the relaxation time for smaller Da remains large, 
while it is quite small for larger Da, as already shown in 
Fig. 3(a). 

Now, we study the mechanism for the observed dras- 
tic enhancement of the relaxation time, by exploring the 
mechanism behind the dependence of the relaxation pro- 
cess on the spatial fluctuations in the case of small Da 
as in Fig. 1. If Ui is distributed almost uniformly, Ai 
is consumed almost uniformly through the entire system 
as in Fig. 1(a). On the other hand, if Ui has large spa- 
tial variations, Ai is consumed locally at the sites with 
large Ui, resulting in an amplification of the spatial fluc- 
tuations, as shown in Fig. 1(d). Then, at sites with 
large Ui, Ai is consumed rapidly so that the reaction 
A + v + 2u ^kBA B + 3u hardly occurs at these sites 
(kAB = 1 >> ksA-)- The decrease of Ui there mainly 
progresses by the reaction u; t ^ Vi which does not con- 
sume Ai. 

On the other hand, with this spatial structure, at sev- 
eral sites that do not belong to the peaks in the pattern, 
Ui and Vi are suppressed. Then, the reaction between 
u and v with the consuming A (whose rate is AiViuf) 
is highly suppressed, even if Ai therein is large. Ai at 
such sites with small Ui is consumed only little by little. 
Hence, a plateau appears in the relaxation of < A >. Af- 
ter the decay of Ui at the site with large Ui is completed, 
Ui at some other site that keeps large Ai starts to be 
amplified by consuming Ai. This process corresponds to 
the re-organization of spatial structure of Ui as described 
in Fig. 1(b) and Fig. 1(d). In this way, several plateaus 
appear successively during the relaxation process. 

Based on the above argument, the drastic increase in 
the relaxation time by the increase of the initial hetero- 
geneity 8 is explained by comparing the speed of the re- 
action consuming Ai with the reaction Ui ^ Vi. Note 
that for a stationary state, Viuf ~ Ui~ Vi is satisfied (by 
disregarding diffusion) if Ai >> 1. When the locally con- 
served quantity Ui + Vi = 2S is large, this stationary so- 
lution satisfies Vi ~ l/ui. Then, the consumption speed 
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of Ai can be approximated by ~ AiViuf <~ Ajitj, accord- 
ing to the estimate for the steady state. On the other 
hand, the order of the speed of the reactions u, ^ Wj is 
given by ~ u,, which is much smaller when < A >3> 1. 
Accordingly, the relaxation process of < A > with the re- 
organizations of m involve fast and slow processes, with 
reaction speeds ~ AiUi and ~ Ui. Thus, the relaxation 
of < A > has several bottlenecks corresponding to the 
slow processes with speeds in the order of ~ m, leading 
to plateaus. In contrast, the relaxation from a state with 
little inhomogeneity involves only a fast process with a 
reaction speed ~ AjUj, since the inhomogeneity in it, is 
not amplified to form a structure. Hence, the relaxation 
time is small when 5 < 5 C . 

Next, we explain the ZJ^-dependence of the relaxation 
time in Fig. 3. As it; increases at a certain site i, the 
resource A diffuses into this site from adjacent sites as Ai 
is consumed at site i. Thus the increase of Ui is further 
accelerated by the use of Aj at the sites adjacent to i. 
Consequently, the peak height of can be much higher 
than is the case when Da = 0. With such increases 
of the spatial amplitudes, the time intervals between re- 
organizations of Ui become larger, resulting in an increase 
of the relaxation time. 

On the other hand, if Da >> 1, the resource A is con- 
sumed faster due to the diffusion of A. In this case, the 
speed of the flow of Ai is higher than that due to the 
reaction u, ^ Vi for sites with large Uj. Therefore, the 
resource A, is consumed continuously by sites with large 
Ui, and a re-organization of the spatial structure Ui no 
longer occurs. In this case, the consumption speed of Ai 
goes up to the level for the relaxation from a homoge- 
neous pattern. Hence, the relaxation time to equilibrium 
for the case Da >> 1 is much smaller than for the case 
with Da < 1, even when S > S c . The peak of the relax- 
ation time in Fig. 3(b) is thus explained. 

Note that the mechanism for the slowing down of re- 
laxation processes proposed here is rather general. Take 
any reaction-diffusion system in which dissipative struc- 
tures are formed by constraining the concentrations of 
some resource chemicals in such a way that their values 
are larger than their equilibrium values. Then consider 
the corresponding closed system, where the dynamics of 
the resource chemical(s) is incorporated. Now the pro- 
posed mechanism for slowing down the relaxation process 
is possible if the following two conditions are fulfilled. 

First, the reaction and diffusion processes of the chem- 
icals that give rise to the non-equilibrium conditions 
should not be too fast compared to those of the other 
chemicals. The diffusion constants as well as the reac- 
tion rates for the consumption of the resource chemi- 
cals should be smaller than those for the others. This 
leads to differences in the time scales of the concentration 
changes and thus, when the pattern formation progresses 
fast enough, the consumption of resource chemicals to 
support the non-equilibrium conditions slows. Then, the 



changes of resource chemical concentrations work as a 
kind of slow variables (or parameters) of the system. 

Second, the consumption of resources should slow 
down due to feed-back from the spatial structure. In the 
present example, the consumption of resources is com- 
pleted soon at sites with higher itj and for other sites 
with very small concentrations of u, the consumption re- 
action progresses only slowly. Hence the overall depletion 
of resources is much slower than in the case with a ho- 
mogeneous concentration of u. 

In general, it is not particularly difficult to satisfy these 
two conditions, and indeed we have confirmed the present 
mechanism by studying some variants of the present 
reaction-diffusion system. 

If only the second condition is satisfied but not the 
first condition, some increase in the relaxation time is 
still observed, but it is neither drastic nor does the re- 
laxation have several plateaus. As an example we con- 
sider the reaction system i)a, + + 2ui a\ + 3«,, 
ii)bi+Ui # bi+Vi, m)ei+Ui # e,-+Cj, iv)f%+Vi # f%+di, 
with diffusion. If we set e, = and c, = c(= const), 
this model is equivalent to a Brusselator0, 0] model 
with the reactions reversed, while it corresponds to the 
Gray-Scott modelUS E3, 

again with the reactions re- 
versed, if bi = and di = d(= const). In this case too, 
Turing patterns are formed if the concentrations of re- 
source and waste chemicals are fixed as a; = a a with 
k « 1 in order to keep the non-equilibrium condition. 
By including the dynamics of ai and by choosing a large 
ai initially, a structure is formed if the initial inhomo- 
geneity 5 is not too small, in the same way as for our 
model above. Indeed, some amplification of the relax- 
ation time is observed here. However, once the initially 
formed structures are destroyed, new large structures are 
not produced. Hence, the enhancement of the relaxation 
time in this case is not as significant as in the previous 
case. 

In this Letter, the relaxation process to equilibrium is 
investigated through a closed coupled chemical reactor 
system. Under certain conditions, we have found that 
the relaxation is drastically hindered once a pattern is 
formed. In addition, we have observed repeated forma- 
tions of patterns. With this itinerancy over patterns, the 
relaxation is further slowed down as compared to the case 
without the structure formation. 

In experimental studies of dissipative structures, the 
system under consideration is usually set to be open in 
order to sustain the non-equilibrium condition. Still, 
even in closed systems, dissipative structures are often 
observed as transients which may last for rather long time 
spans (recall for example the Belousov-Zhabotinsky re- 
action in a petri dish). By choosing a suitable reaction 
system, it will be possible to demonstrate the present 
enhancement of the relaxation time due to the transient 
dissipative structure. In relationship with the possibil- 
ity of an experimental observation, it should be noted 
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that the slow relaxation process we studied here does 
not progress gradually, but repeats residence at several 
quasi-stationary states corresponding to certain struc- 
tures. The quasi-stationary states associated with the 
successive spatial structures act as bottlenecks for the 
relaxation process. 

In complex reaction systems, with more chemical com- 
ponents, the relaxation process could further be slowed 
down. For example, assume that A and B in our model 
are synthesized by lower-level resources A 1 and B' , and 
that these reactions also satisfy the mechanism demon- 
strated in the present paper. By a hierarchy of such 
reactions, the relaxation time is then expected to be fur- 
ther increased. This may provide some insight into why 
a cell system can maintain a non-equilibrium state over 
a huge time span. 
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for Scientific Research from the Ministry of Education, 
Science and Culture of Japan (11CE2006). 
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FIG. 1: Typical temporal evolutions of it; (left), A t (right) 
for D v = 250, Da = 0, and S = 4. (a)S = 0.1 and (b)S = 4.0 
plotted until < A > becomes smaller than 0.001Ai ni . (c) The 
time evolutions of < A > corresponding to (a) and (b). (d) 
Three typical snapshots of spatial patterns in (b) which are 
plotted at the time step shown by the arrows. 




-0.4 



-0.6 I 

100 
F 



FIG. 2: Time course of (< A >' (t),F(t)) obtained from the 
same simulation of Fig. 1(b). See text for the definition of 
< A >' and F. 
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FIG. 3: (a) The average relaxation time T, plotted as a func- 
tion of 8, for Da = 0, = 0.25, and = 25, and (b) T as a 
function of D A for 8 = 0.1 and = 4. D„ = 250. 



